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SUMMARY 


This investigation is intended to improve the numerical practicality of a matched 
asymptotic expansion approach for the computation of unsteady three-dimensional 
airloads on a helicopter rotor. The original method as suggested by Van Holten has 
previously been evaluated and proven to be a comprehensive and accurate analysis for 
flight conditions conducive to linear flow phenomena. This effort to decrease the 
computational requirements of the original analysis utilizes a discretized representation 
of the doublet strength distribution and helical streamlines. The continuous variation of 
the doublet strength has been approximated by piecewise constant or piecewise 
quadratic distributions, and the helical trajectory of a fluid particle has been 
approximated by connected straight line segments. As a direct result of these 
simplified representations the computational time required for the execution of a 
typical flight condition has been reduced by an order of magnitude with respect to the 
requirements of the original analysis. Airloads which have been computed using the 
discretized method for a two-bladed model rotor and a full-scale four-bladed rotor are 
in close agreement with measured results and airloads from the original asymptotic 
analysis. For conditions characterized by significant rotor/wake interaction the 
piecewise constant representation requires a reduced azimuth spacing to maintain 
acceptable accuracy. 


INTRODUCTION 

The problem of estimating airloads on helicopter rotor blades can be approached 
by a variety of approximate methods. One such approach, put forward by Van Holten 
(refs. 1-4), uses an acceleration potential description of the flow field and a matched 
asymptotic expansion technique. Under the assumption of incompressible potential flow 
the unsteady three-dimensional airloads on a rotor blade in forward flight are calculated 
to a consistent order of approximation in terms of the aspect ratio. 

A study has been conducted (ref. 5) to examine the theoretical basis and 
computational feasibility of the method, and to evaluate its performance and range of 
validity by numerical comparison with experiment and other approximate methods. The 
study concluded that, within the restrictions of linear theory (i.e., small disturbances), 
the Van Holten approach does lead to a valid description of the rotor flow field and a 
systematic determination of the airloads on the rotor blade. It was also found for flight 
conditions involving significant blade/wake interaction effects, the agreement between 
computation and measurement is poorer than in other cases. 

The analysis in Van Holten's approach leads to an integral equation for the blade 
doublet strength distribution (eq. (8), ref. 5). This is solved using a collocation 
technique, which consists of assuming the unknown function to be made up of a 
combination of suitable spanwise modes and azimuthal harmonics, and satisfying the 
integral equation at an equal number of points distributed over the rotor disk. The 
result is a set of linear, simultaneous algebraic equations. However, setting up the 
equation at any collocation point requires integration with respect to azimuth of the 
individual assumed mode combinations. This numerical integration makes up the bulk of 
the total computation required for the solution. Under conditions of low forward speed 
and low inflow (when a larger azimuth range must be covered with a finer integration 
step) and/or a larger number of blades, the computation time is significantly increased. 
As an example, for a two-bladed rotor at an advance ratio of 0.3, the computer program 


takes about 250 seconds to execute on a CDC 6400 computer. For other conditions, the 
time would vary approximately in direct proportion to the number of blades and in 
inverse proportion to the forward speed (advance ratio). In addition, when interaction 
effects are judged to be significant (possible even at moderate to high forward speeds), 
a smaller integration step will have to be used. 

It would therefore be desirable to seek a simplified computational scheme for the 
basic asymptotic approach that would lead to significantly lower computation time 
without sacrificing any of the essential features of the asymptotic method. One 
possibility is to consider the analogous situation in the vortex representation of rotor 
wakes. In the vortex approach, the continuously varying bound circulation on the blade 
generates a wake in the form of a helical vortex sheet of continuously distributed 
trailing and shed vorticity. Calculation of the induced velocity due to such a wake 
would require double integration over the wake surface, which can be time-consuming. 
In practice, this problem is often overcome by assuming that the variation of the blade 
bound circulation over the span and the azimuth takes place in discrete, finite steps. 
This results in a wake of discrete trailing and shed vortex elements. Since the velocity 
induced by a straight vortex element can be analytically expressed, the induced velocity 
due to the entire wake can be written as a summation of analytical expressions 
representing the contribution of individual trailing and shed vortex elements. This 
results in considerable reduction of computational time over the exact numerical 
integration. 

It appears that a similar simplification could be achieved by approximating the 
continuous variation of the doublet strength distribution in the asymptotic approach. It 
is the purpose of the current study to develop a discretized representation for the 
asymptotic method and compare its performance with the original computational 
scheme and its results with measured data for the flight conditions considered in 
reference 5. This report describes the details of the simplified scheme and discusses 
the computational results. Detailed analytical expressions are presented as Appendices. 

SYMBOLS 


A 

a 

, o 

a r b i 

B 

b,c 

c lj’ c 2j’ c 3j 

d 

d 0’ d r d 2 

d 


d n0 ,d nl ,d n2 

F 2’ F 3 


aspect ratio 
coning angle 

blade flapping coefficients 
number of blades 

semi-chord and chord, repsectively 

coefficient of piecewise quadratic representation (eq. (Al) ) 

distance between fluid particle and collocation point 

2 

coefficients of quadratic expression for d in terms of 0 

distance between fluid particle and point on the chord at the same 
spanwise location 

2 

coefficients of quadratic expression for d n in terms of 0 

terms appearing in the regular part of the near field solution (eq. (E 5) ) 


2 


G°,G C ,G S 
J jn’ ]n 


g; 


I C 

1 


harmonic coefficients in the Fourier expansion of e in terms of HV 
(eq. (ED) b 

function representing the doublet strength distribution along the blade, 
and the basic unknown in the problem 

value of g at the beginning of the j n spanwise segment 
i th integral expression for the common part, Appendix B 
i ** 1 integral expression for the far field, Appendix B 
i th integral expression for the near field, Appendix B 


L 

l 

M 

m 

P 

R 0 ’ R i 

R i’ R s 

r 

r b 

r bo 

r o ,r r r 2 


s 

s. 

J 



U 


T 


v. 

10 


w 


Aw 


x,y,z 

Wb 


total lift on one blade (eq. (E 10) ) 

lift per unit span, sectional lift, (eq. (E 7) ) 

moment of lift distribtuion about the rotor hub (eq. (E 11) ) 
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(eq. (E 8 ) ) 

perturbation pressure 

root and tip radius of blade, respectively 
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spanwise location of fluid particle 
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rotor coordinate systems (fig. 1 ) 
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chordwise location of center of pressure (eq. (E 9) ) 

coefficients of linear expression for in terms of 0 (eq. (B 2) ) 

coefficients of linear expression for in terms of 0 (eq. (B 2) ) 

blade inertia coefficient for flapping (Lock number) 

linear blade twist, root pitch angle-tip pitch angle 

azimuth position with respect to reference blade, ( 4^ - 4^ ) 

collective pitch angle at blade root 

nondimensional rotor inflow, ua + v. 

r 10 

nondimensional rotor forward speed 

R 0 ,R 1 
air density 

plane elliptic coordinates 

meridional angle in cylindrical coordinates 

azimuth position with respect to downwind reference line 

azimuth position of collocation point 

azimuth separation of j blade from reference blade 

rotor angular velocity 

DISCRETIZED ASYMPTOTIC REPRESENTATION 


The essential features of the asymptotic approach are retained in the discretized 
representation, i.e., the blade pressure field is obtained in composite form as the sum of 
far field, near field and common part expressions in such a way that the composite field 
reduces to the near and far field solutions in the near and far field regions, respectively. 
In view of the approximation^ being considered, it is not known if it would be consistent 
to retain the terms of 0(1/A ) that are present in Van Holten's analysis (ref. 1). For the 
sake of simplicity in computation these terms are dropped in the present study, so that 
the discretized representation is of 0(1 /A), comparable to standard lifting-line 
formulations. 

The approximations to be introduced are: 

(1) The continuous variation of the unknown blade doublet strength distribution 
over the range of blade span and azimuth is replaced by simple, piecewise continuous 
functions over suitable subintervals of span and azimuth. 

(2) The helical trajectory of a freestream fluid particle relative to the blade is 
replaced by a series of connected straight line segments. 

Over the azimuth, the doublet strength variation is assumed to be piecewise constant, 
i.e., its value is constant over each subinterval of azimuth. Over the span, it is 
approximated by a piecewise constant or a piecewise quadratic variation, the latter 
providing a more accurate spanwise representation at the cost of more complicated 
algebra. Appendix A describes the manner in which these representations are 
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formulated, while typical representations are illustrated in figures 2(a) and 2(b). 

The basic problem is the calculation of the vertical velocity induced on the blade 
by its pressure field, as given by the relation 



- 00 


( 1 ) 


The rotor coordinate systems are shown in figure 1 and ^ is the azimuth position of 
the reference blade. The blade pressure field is written in composite form as 


p = p, + p - p 

K K far *near ^common 

By appropriate construction and matching (ref. 1), all three components can be 
expressed in terms of the doublet strength distribution along the lifting line. 

In Van Holten's analysis, the far pressure field is the field of a dipole line along 
the quarter-chord location and is expressed as a series of associated Legendre functions 
in terms of prolate spheroidal coordinates. The near pressure field is the local two- 
dimensional field of the section and is expressed in plane elliptic coordinates. The 
common part corresponds to the behavior of the far field solution in the near field or, 
equivalently, that of the near field solution in the far field. With the form of the 
pressure field established, the induced velocity is calculated by numerical integration of 
the composite pressure gradient, for which purpose the above expressions are 
convenient. In fact, one of the merits of Van Holten's approach is that the composite 
pressure field is given by a direct expression free of integrals, so that the induced 
velocity calculation requires only one numerical integration. 

However, these expressions are not convenient for the purpose of applying the 
proposed approximations. To 0(1 /A), the far field is the field of a dipole line along the 
blade midchord location, and can be written as the spanwise integral of a distribution of 
three-dimensional doublets. If the approximations described above are introduced, it is 
found that the calculation of the induced velocity due to the far field reduces to a 
double summation (over spanwise and azimuthal segments) of fairly simple integrals 
that can be analytically evaluated. Likewise, the near field is written as the integral, 
over the chord, of a two-dimensional doublet distribution, of strength proportional to 
the surface pressure differential. The chord is divided into a suitable number of 
segments over each of which the pressure differential is assumed constant, at its 
average value on that segment. With this simplification, the induced velocity due to the 
near field also becomes a double summation (over chordwise and azimuthal segments) of 
analytical expressions. The contribution of the common part, which is the field of a 
single two-dimensional dipole, presents no problems. The complete details of the 
various steps pertaining to the above calculations are presented in Appendices B and C. 
Appendix B describes the result of approximating the trajectory and also presents a list 
of integral expressions to be used for the induced velocity calculation. Appendix C 
derives the summation expressions for the induced velocity due to the far field, common 
part, and near field, for both the piecewise constant and the piecewise quadratic 
representation of the doublet strength variation along the span. Appendix C also 
demonstrates that the expressions derived for the far field and the common part 
correctly cancel in the vicinity of the blade, thus verifying their asymptotic character. 

In summary, the effect of the discretized representation on the problem is to 
reduce the induced velocity integration of equation (1) to a summation of analytical 
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expressions. The normal velocity boundary condition is then applied by setting this 
induced velocity equal to the normal velocity on the blade surface due to blade motion. 
The form of this equation is given in Appendix E (eq. (E 2) ). The unknowns to be solved 
for are the discrete values of the doublet strength at the spanwise segment midpoints, 
at discrete azimuth locations. However, for final presentation as output, the discrete 
variation with azimuth would have to be fitted by an interpolation curve. This is 
accomplished by substituting into equation (E 2), for the azimuth variation of each 
midpoint value, a harmonic interpolation formula given by (E 3). It must be noted that 
this does not imply a continuous azimuthal variation in the induced velocity calculation, 
since the doublet strength is still held constant over each azimuth segment. The 
advantage of the substitution is that the interpolation coefficients are obtained directly 
and available for calculation of output data at desired azimuth locations. Expressions 
for the various quantities calculated for presentation as output are given in Appendix E. 

Examination of equation (E 2) shows that the near field induced velocity contains 
a contribution from the regular solution, proportional to the blade motion parameters. 
If these parameters are considered known, this contribution should be taken together 
with the blade normal velocity term. However, if any or all of the blade motion 
parameters (collective pitch, coning angle, cyclic pitch coefficients) are considered 
unknown, additional equations must be generated to solve for them. This can be done by 
using the following conditions: 

(1) Azimuth average of the total lift due to all blades should equal the known 
rotor thrust. 

(2) Flapping moment equilibrium should exist about the rotor hub. 

If only first harmonic flapping is considered, the second condition yields three equations 
(zeroth harmonic, first harmonic cosine and first harmonic sine components). These 
additional equations are also listed in Appendix E, for the piecewise constant and 
piecewise quadratic representations. 

DISCUSSION OF RESULTS 

The discretized representations which have been incorporated in the analysis have 
been numerically evaluated to ascertain their potential accuracy. All analytical details 
of this evaluation are presented in Appendix D. Figure 2 illustrates the piecewise 
constant and quadratic representations of the spanwise distribution of doublet strength 
for a typical condition. 

For the near field chordwise approximation, an airfoil with a steady two- 
dimensional pressure distribution is considered. The induced vleocity of a fluid particle 
travelling parallel to the chord is calculated for various vertical distances from the 
airfoil. The calculation is performed using both the exact solution and the 
approximation. The results are tabulated in Appendix D (table Dl) and plotted in figure 
3. It can be seen that the approximations with three segments and five segments 
across the chord are acceptably close to the exact solution, even at very small vertical 
distances. This simple example shows that, as far as the induced velocity calculation is 
concerned, even a relatively crude representation of the surface pressure distribution is 
sufficient. 

For the spanwise approximations of the far field solution, a finite wing is 
considered with a spanwise distribution of a form which is typical of a rotor blade 
distribution. The induced velocity due to this distribution is calculated for a fluid 
particle travelling parallel to the chord, at various spanwise locations and vertical 
distances. The calculation is performed using the exact solution as well as the 
piecewise constant and piecewise quadratic approximations with three and five 
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spanwise segments. The results are tabulated in Appendix D (tables D2 and D3) and 
plotted in figure 4. With the three segment model, it can be seen that both 
approximations deviate from the exact solution when close to the wing, the deviations 
being more marked near the loaded tip and generally greater for the piecewise constant 
representation. With five segments, the piecewise quadratic representation is nearly 
identical to the exact solution while the piecewise constant results still show significant 
deviation very close to the wing. It may be noted that the exact and approximate 
results tend to merge with increasing distance from the wing, as is to be expected. It 
may also be noted that the results for different spanwise locations tend to merge with 
increasing vertical distance. 

Airload computations have been carried out for the same experimental cases 
considered in reference 5, viz, (1) a two-bladed teetering model rotor at forward 
speed/tip speed ratios (y ) of 0.08, 0.15, 0.29 (ref. 6), (2) a four-bladed articulated full- 
scale rotor tested in flight, at y =0.06, 0.13, 0.29 (ref. 7), and (3) the same four-bladed 
rotor tested in a wind tunnel at y = 0.29, 0.39, 0.45 (ref. 8). The geometric and flight 
conditions are listed in table 1, which is reproduced from reference 5. The measured 
results are compared with computations using the piecewise constant (p.c.) and 
piecewise quadratic (p.q.) representations, as well as Van Holten's computational 
scheme. Before proceeding with a discussion of these comparisons, some comments 
regarding the computations are in order. In the original scheme, the numerical 
integration is carried out with a 5-point Gauss-Chebyschev rule over a suitable 
subinterval of azimuth, to be properly chosen for accurate computation. In the 


discretized representation, there is no numerical integration but a choice has to be 
made with regard to a suitable azimuth subinterval over which the trajectory is 
straightened. In both cases, the smaller the azimuth subinterval chosen, the more exact 
are the computations. In order to keep the computation economical and at the same 
time achieve some of the accuracy of small azimuth spacing, the procedure adopted is 
to use a "normal" spacing whenever the fluid particle is not close to the blade and a 
"reduced" spacing whenever it is close to a blade. Reference to figure 1 shows that the 
trajectory locations at which the fluid particle is directly over a blade are 


characterized by = 0. For each collocation point (r, , V, ) these azimuth positions 
are determined in advance using equation (B 1). Durmg tne azimuth integration (or 


summation), a reduced spacing is used in the vicinity of these locations. In the 


computations carried out here (unless otherwise mentioned) the normal and reduced 
azimuth intervals used were 15 and 5° for the Van Holten scheme and 30° and 10° for 


the discretized representation. Although results of computing with the original scheme 
were reported in reference 5, these results were recomputed for the present study using 
normal and reduced spacing as above, and with spanwise collocation points at 
r/Rj = 0.30, 0.55, 0.75, 0.85, 0.95. Due to these changes, some differences will be noted 
between the curves presented here and the corresponding ones in reference 5, but the 
differences are generally small with one exception which will be pointed out later on. 
For the discretized representation five spanwise segments were used, with end points at 
r/R, = RJR., 0.5, 0.7, 0.8, 0.9, 1.0, and five chordwise segments with end points at 
x/b = - 1.x), -0.9, - 0.6, 0.0, 0.5, 1.0. 


Results for the variation of total blade lift as a function of azimuth position for 


Case 1 are shown in figure 5. It can be seen that for all three forward speeds, the 
results of the discretized representation are quite close to the original scheme, with the 
p.q. representation being generally a little closer than the p.c. representation. 

For Case 2, the results are illustrated in figure 6 and, as may be anticipated from 
the computations reported in reference 5, there is greater variation here, mainly for 
y = 0.13 and y = 0.29. At y = 0.13, the original scheme results in a considerably wavy 
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TABLE I. - GEOMETRIC AND FLIGHT CONDITIONS FOR THE 
EXPERIMENTAL CASES 
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distribution which is not seen in the results of the discretized scheme, especially with 
the p.q. representation which remains close to the measured curve except near the 
advancing blade position. At p = 0.29, there is once again a deviation near the 
advancing blade position, with the p.c. results being particularly bad in this region. The 
p.q. results are generally close to the original ones. However, both results from the 
discretized scheme tend to overestimate the lift in the disk trailing edge region (around 
Y b = 0°). 

The results for Case 3 are presented in figure 7. It was observed in reference 5 
that the original results showed a tendency to overestimate the lift near the downwind 
azimuth position and correspondingly underestimate it near the upwind position. This 
tendency is also present in the results from the discretized representations. This may 
be due to the increasingly important effect of radial flow (along the blade span) in these 
regions at moderate to high forward speeds. The influence of radial flow is primarily on 
the spanwise development of the blade boundary layer, increasing with forward speed, 
and must be accounted for empirically. Otherwise, the p.q. results are acceptable and 
show the same trend as the measured curves. However, the p.c. results show 
particularly significant deviations for p = 0.29 and p = 0.39, near the advancing blade 
position and in the disk trailing edge region. It was noted earlier that the changes made 
in the original scheme produced only small variations with one exception. This 
exception is p = 0.39 of Case 3. Comparison of figure 19 of reference 5 with figure 7(b) 
of the present study reveals that the latter variation is much better, particularly in the 
absence of the large peak near the disk trailing edge that is present in the former. This 
difference is surprising, at first glance, because the changes made do not seem that 
important. Apart from using a constant azimuth spacing, the original version of the 
computational scheme also used spanwise collocation points at r/Rj = 0.40, 0.55, 0.75, 
0.85, 0.95, differing from the present version only in the innermost point. Both 
programs use the zero lift condition instead of the normal velocity boundary condition 
at a collocation point whenever the local effective freestream speed U^.< 0.1. 
However, if r/R^ = 0.4 is used, there is a collocation point near the retreating 1 blade 
position that has 1 = 0.105, which is small yet large enough to escape the zero lift 
condition. It is apparently this particular point that is the source of the trouble, for 
when it is replaced by another collocation point with U~> 0.2, the resulting distribution 
is quite regular (like figure 7(b) of the present study). This situation does point out the 
need to use care in the choice of collocation points (avoiding those points with small 
positive values of U^.) and use of close azimuth spacing in those cases where 
interference effects are significant. 

Spanwise distributions of sectional lift, at various azimuth positions, are plotted in 
figure 8 for p = 0.29 of Case 1, and in figure 9 for p = 0.29 of Case 2. General 
comments on these curves are much the same as those made for the original scheme in 
reference 5. Agreement is acceptable as long as the measured curves do not show any 
sharp variations, as in the case of close interaction with a tip vortex. When sharp 
variations do occur, they are not evident in the computed curves. There is also general 
deviation from the measured curves in the tip region. In addition, the failoff to zero at 
the tip is better with the piecewise quadratic represenation, presumably because it 
provides a better approximation to the actual curve than the piecewise constant 
representation. 

In discussing the total lift variations it was noticed that the p.c. representation 
led to particularly significant deviations for three flight conditions , viz, p = 0.29 of 
Case 2, p = 0.29 and p = 0.39 of Case 3. To test the possibility of improvement with 
smaller azimuth spacing, these three cases were computed again with a constant 
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azimuth subinterval of 5 and the results are shown plotted in figures 10 to 12. It can 
be seen that there is considerable improvement in the results for the p.c. 
representation, while the p.q. results are only slightly changed. This is a clear 
indication that, for flight conditions involving significant interaction effects, the results 
from the p.c. representation are sensitive to the azimuth spacing used and should be 
calculated with a small spacing. 

Since the primary objective of this investigation is to improve the numerical 
practicality of the original asymptotic execution time requirements, it is essential that 
these requirements be examined. Table 2 presents results suitable for comparison of 
the different representations for two typical rotor configurations. In addition to the 
original continuous representation suggested by Van Holten and the piecewise constant 
and quadratic representations of this study, there is also presented data for the 
segmented lifting line model with a discrete vortex wake. This was the only linear 
method considered in reference 5 which provided airloads of an accuracy which could be 
compared with the asymptotic method. It can be observed from table 2 that the 
piecewise constant representation reduces the execution time requirement of the 
original continuous representation by a factor of seven. A reduction to almost one fifth 
is attained for the piecewise quadratic simulation. As indicated in the table, these 
figures are based on computational azimuth intervals of 30° normal and 10° reduced 
spacing. In cases where much shorter intervals are required for the piecewise constant 
representation the computational requirements are proportionately increased as 
illustrated in Table 2. It should be noted that the computational efficiency of the 
piecewise constant representation is equivalent to the less comprehensive segmented 
lifting line with a discrete vortex wake when the azimuth intervals are the same. 


TABLE 2. - COMPARATIVE COMPUTER EXECUTION 
TIMES FOR TYPICAL CONDITIONS 


Computational 
azimuth intervals, 
normal/reduced 

Method of 
analysis 

CDC 6400 
execution time, 
sec 

Case 1 y = 0.29 

1 5°/ 5° 

Original continuous 

231 

30°/ 10° 

Discrete vortex wake 

30 


Piecewise constant 

33 


Piecewise quadratic 

49 

Case 2 y = 0*29 

1 5°/5° 

Original continuous 

467 

30°/10° 

Discrete vortex wake 

60 


Piecewise constant 

65 


Piecewise quadratic 

100 

o 

o 

Discrete vortex wake 

265 


Piecewise constant 

265 


Piecewise quadratic 

415 
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CONCLUDING REMARKS 


The asymptotic approach developed by Van Holten is a suitable method for rotor 
airload calculation, within the scope of linear theory. However, in spite of the fact that 
only a single numerical integration is required to calculate the induced velocity, 
significant computation times may be required under certain conditions. It is possible 
to reduce the computation time by making two approximations, viz, replacing the 
continuous variation of the doublet strength distribtuion along the blade span by a 
piecewise continuous variation, and replacing the continuous helical trajectory of a fluid 
particle by a succession of connected straight line segments. 

In the present study, such a discretized representation has been developed for the 
asymptotic method, using either a piecewise constant or piecewise quadratic variation 
of the doublet strength along the span. Computations have been carried out for the 
case of a two-bladed, teetering model rotor and a four-bladed, articulated full-scale 
rotor, and the results compared with both measurement and the original computational 
scheme. In general, when interaction effects are not significant, the simplified scheme 
agrees well with the original results, with the piecewise quadratic representation being 
slightly better. When interaction effects are significant, the piecewise constant scheme 
yields poor results but is found to improve upon using smaller azimuth spacing, while the 
piecewise quadratic scheme continues to compare well with the original results. 
Computationally, the discretized representation shows considerable improvement over 
the original scheme; under conditions where interaction effects are not significant, the 
piecewise constant scheme requires only about one-seventh and the piecewise quadratic 
scheme about one-fifth of the computation time required for the original scheme. 
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APPENDIX A 


PIECEWISE CONTINUOUS REPRESENTATIONS 


The continuous variation of the dipole strength function, g, over the blade span is 
approximated by dividing the span into 3 segments and replacing the actual curve by a 
series of simpler curves, continuous over each segment but discontinuous at the segment 
boundaries. The points marking the spanwise division are labeled s., j = 1,2,. . .(3+1) 
with Sj being the root and s^ + j the tip of the blade. Points r., j = 1,2. ... 3 denote 
the midpoints of the segments, where the boundary conditioh of normal velocity is 
applied. The length of a segment is As-, j = 1,2. . .3 and this is allowed to vary along 
the span, generally being chosen smallef near the tip to achieve a better representation 
of the rapid variation of blade loading. 

The simplest representation is a constant value over each segment and this is 
shown in figure 2. The spanwise variation can be represented much better with a 
second-degree curve over each segment as illustrated in figure 2. However, setting up 
such a piecewise quadratic representation involves more complicated algebra, the 
details of which are given below. 

The quadratic curve is constructed to satisfy the following requirements: (1) at 
the midpoint, r., it must have the actual functional value, g(r.), (2) at its end points, s. 
and s. + j, it mus\: equal the values on the adjacent segments, Oj at its end points, it musl! 
also 4qual the slopes of the adjacent segments. It must be noted that, for the segments 
at the ends of the blade, the curve is allowed to go to zero at the root and the tip, 
without any constraint on the slope at these points. Since each curve requires 3 
constants for its definition, a total of 3 3 constants must be determined from the 
available conditions, which are listed below. 

Number of equations for midpoint values = 3 

Number of equations for end point values = 3+1 

(including zero values at root and tip) 

Number of equations for end point slopes = 3-1 

(excluding root and tip) 


This adds up to a total of 3 3 equations, so that it is possible to set up a piecewise 
quadratic approximation which is completely determinate. 

The curves will have to be determined in terms of the midpoint values and 
segment locations. To begin 
values assumed known. If C 
conditions to be satisfied are 

2 

C, . + C-. s. + Ci- s. = g. 
lj 2 J J 3j j 5 j 

C lj + C 2j s j+l + C 3j s j+l = g j+l 

which can be solved to give 


C, . + C-. r. + Co. r. = g(r 
lj 2) j 3) ) 6 


with, the slope conditions are ignored and ||ie end point 
j., C 2 j, C^. denote the constants for the j curve, the 
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c ii = 2 [ r j 5 j ® j+ 1 * r j Vi g f 2s i Vi g(r i > ] A s f 

c 2j = - 2 [ (r j * V «j t i * <r j + Vi 1 V " r i g<r j^ is j 2 (A “ 
C 3j = 2 [s j+ 1 + Sj - 2 8(rj>]/ as f 


Now the slope conditions can be used to elimate the end point values 
j = 1,2. . . (3+1). At any end point, say Sj + j, the equality of slopes leads to 



+ 2 C,. s. , 
3j j+1 


2 C, . , s. , 
3,j+l j+1 


which, on substitution, becomes 

U/As.) gj . 3(l/4 Sj ♦ lMs j+1 )g jtl ♦ UMs jtl )g jl2 = ♦[g(r j )/45j + g(r j+ ,)Ms j+ ,] 

for j = 1,2. . .(3-1) (A 2) 

This tridiagonal system of (3-1) equations can be solved by the following recursive 
scheme, which consists basically of the forward elimination and back substitution of the 
Gaussian elimination method. 


3 2 = 3(l/ASj + 1 /As 2 ) 

Y 2 = 4 [§ (r p/ASj + g(r 2 )/A sj 

3j = 3(1/As._ 1 + l/ASj_ 2 ) - 1/As?_|3j_ 1 

Vj = {*D str i-i >/As j-i + « (r i )Ms j] 

- Tj_ 1 /A s j_ 1 1 / Bj (j = 3, k. 

g 3 = Yj 



gj = Yj ~ gj + l/Sj As j (j = 3-1, 3-2. . .2) 
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This determines the (3-1) end point values in terms of the 3 midpoint values, so that the 
quadratic components C.., C--, Cy are completely defined. Within a segment, the 
function can be written as^ ’ ‘ 


2 

g = C, . + C 0 . s + C,. s 
6 lj 2j 3) 

= 2 (s - s. ,)(s-r.)g. + (s - s.) (s - r.) g. . 

1_ j+1 j °j j j j+1 

- 2(s - Sj) (s - Sj + |) g(r)J /A (A 4) 

Calculation of the total blade lift and the moment about the hub requires the 
spanwise integral of g and its moment about the hub. These integrals are given below. 

Piecewise constant representation: 


i 

/ 


g(s) ds = 


E as. 
j=l > > 


I 

I 


g(s) s ds = 


Z g( r j) fjA s- 
j=l 111 


Piecewise quadratic representation: 


(A 5) 


/ 

S 

1 

/ 


g(s) ds = 


g(s) s ds 



s. 

J 


6 


(A 6) 
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APPENDIX B 


TRAJECTORY APPROXIMATION AND 
LIST OF INTEGRALS 


The helical trajectory of a freestream fluid particle, ending at a^ollocation point 
(r b ^ R L’ ^bo^ on a re * erence blade, relative to the axes fixed to the j n blade, is given 
by ?he iollowing equations. 


x b /Rj = (r bQ /R^sin(8 + A Y) + y0 sin ( 0 + Y bo + AYJ 

y b /R 1 = X0 


(B 1) 


r b /R l = (r bo /R l )cos(0 


AY.) + p0 cos (0 + Y fao + A Y.) 


where 0 = (Y b - Y b ) is the azir^juth relative to the reference blade position and AY- is 
the azimuth separation of the j n blade from the first (reference) blade. For uniforrhly 
separated blades, 


AY. = 2 tt ( j - 1)/B 

The approximation used here divides the continuous helical trajectory into a 
succession of straight-line segments, each connecting the initial and final points of the 
helical path over a subinterval of azimuth. The coordinates of each of these trajectory 
segments can be written as 


V R i = x i + x s 0 

V R 1 = Y i + Y s 6 


r b / R , = Rj . R $ 8 


(B 2) 


If 0j and 6 2 are the ends of the azimuth interval, the intercepts and slopes are given by 

X i = (x bl - R l X s 9 1 ) ' R 1 
X s = ^ x b2 " x bl^ / R i ( 0 2 " e p 


Y. = 0 

l 

Y s = * 


R i " (r bl ' R l R s 6 1 ) 1 R 1 


R s = ^ r b2 " r bl^ R i^ e 2 * 


(B 3) 
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The integrals along the fluid particle trajectories which appear in Appendix C contain 
the following expressions which are written in terms of the straight-line notation as 


2 2 2 / - 0 a 2 Xri 2 
r = x b + y b = (r o + 2r l 6 + r 2 9 )R 1 


where 


r = XT 


r, = 


r-i = 


X. X 

1 s 

2 2 
X + 

s s 


Letting u = 


<s - r b ) 


d 2 = r 2 ♦ u 2 = (d ♦ 2d,0 + d,0 2 ) R? 

O 1 l 1 


where 


d = X 2 + (s - R.) 2 
o 11 


= 


d, = X. X - (s - R.) R 
1 is is 

2 2 2 
d 2 = X 2 * X 2 * R 2 

(x-x b )2 * y 2 = (d no * 2d n ,0 * d n2 0 2 )Rf 


where 


nO 


(x - X.)' 


nl 


n2 


- (x - X ( ) X $ 

X 2 + Y 2 
s s 


The primary advantage which has been achieved by introduction of the straight- 
line segment approximation for the particle trajectories is evident in the above 
"distance" expressions. It can be noted that they are all quadratic expressions in 0. As 
a direct result the integrals of the pressure gradient which are described in Appendix C 
can be analytically evaluated. These . integral expressions are listed below as 
determined from reference 9. The symbolic parameter, t, is defined as 


t = dj + d20 + dVcl^" 
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Far field integrals: 


. f f d0 £n t 

'i : J d : r 


f ede 

J d = 


d - d iV 


i f _ f d e f de 
*3 ~ J d(d - u) " J d(d + u) 




t + (R d - uVcTZ) 
s i 


x.x 


- tan 


t - (R g d - uVcTpj 


,jl^ 


I * - j * 

S " l 3 


y b = y o + X6 


f 9 d 9 
" J d(d + u) 




y b = y o + A6 


y •* 0 
3 o 


2X(R d-uVdT) X I, f 2y.u 


V " u y 2 3 

2 ' X.X 


X.X . . 2 

l tdr 


2X s (R s d-u^)(X s x bt Xy b ) 
Xj X 
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T f 9 _V 

*8 = 3 y r 


y -*• 0 
} o 


2 R $ ad+ y b Vd^) 

tdVd- 


2 yK u 


dr 


, f f . f 
+ A - r ^ ly 


I, 


f _ f e 2 de f e 2 d e 
9 ~ J d(d - u) ‘ J d(dT 


ST 


[ 2 (s - Rj) 1 / - 2R S I 2 f - r o I 3 f - Zrj l/ ] / 


, * - ( 

10 J , a 


d(d - u r 


f— 

^ d(d + 


de + f de 


u)‘ 


1 


r 2 tX . 2 X 2 VdT 

i 2 


+ t(2u + r" 


-J 2(d 1 2 - d 0 d 2 ) (dVd^ - R g u) 
I 3 ,Va ?[ r l R s * r 2 <s - R i ) ] } 


t f f Q d9 

11 " J d(d - u > 2 * 


r 9 d 
J d(d + 


u)' 


= - (2d/r 2 + r A 0 + R g 0 / r 2 
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f r e 2 de r e 2 d e 

12 J d(d - u) 2 J d(d + u) 2 


= ^21^ + 2{s - R.) I 3 f - 2R s I, 


" r 0 *10 


_2r i 'ii*]/ 


f _ f 9 3 d 9 + f e 3 d 9 

J d( d + 


, f r e de 

13 = Jd(d-u) 2 


= [2I 2 f + 2(s - Rj) I 3 f - 2R s I 9 f 


" r o *11 ~ 2r 


1 * 12 *] 1 


f 9 d e 
J d(d + u) 


£,n[2d,t (d + u)]] 


2 VT 2 s' 


2r i . -1 { * - < R s d - uvap 
r 2 l x i x l I I Xj X I 


19 



= 2 h l + (s - R i> >3 f - R s 1 5 

i f _ r ode f 9d_e 
17 ~ J d-u + J d+u 

= 2I 2 f ♦ (s - Rj ) I, f - R I, f 

hi = / [ y b fi ' n(d + u) ] d 9 

= 9 £n(d+ u) + R I- f - r. I* - X 2 I* 
s Z 1 s 15 
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Some of the above integrals are not valid for the trajectory segment which ends at 
the collocation point for which X. = 0. For this circumstance the following integrals are 
used. 
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Common part integrals 



= (r j + r 2 0 + r 2 r 2 1^) / 2X. 2 X 2 r 2 

, c r ede 

3 = J S 

= - (1 + 2r 1 r 2 I 2 C ) / 2r 2 r 2 
c f e 2 de 

■4 _ / 4 


= (r o r2l 2 C - 0)/r 2 r2 

c _ f e 3 de 

5 = J r 4 


= (5,n r - r . I. C - r I, c 
11 o z 3 


- 2r ] r 2 ^ C) / r 


2 

2 
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/ tl 

9 de 

= (e 3 - 3r o r 2 I^ c - r 2 I 5 C )/r 2 r 2 
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= (£nd -d . I, n )/d 
n nl 1 r 


n _ r e 2 d-e 
3 J d 2 

n- 


<« - d no [ l" - 2d nl l 2 n > 1 d 


n2 


I 


e 3 d 9 

4 ■ 1 d 2 

n 


"■/ 


- <« 2 - 2d no l 2 - 4d nl 1 2d n2 


APPENDIX C 


INTEGRATION OF PRESSURE FIELD 


The composite pressure field of the blade, to 0(1/A), can be written from 
equations (6) and (7) of reference 5 as 


2„2 

pft R, 


r p di P i + 

1(1 - 0 g Sin Y , 

[ g sin d) 

Lpa^jJ 

L 2A pft 2 R 2 r/R l . 

L p s} 2 Rj cosh n * cos * 


+ (F, + (r b /R | ) F^e~^ sin <}> 


(l-C) 


for points within the blade span, and for points outside the span 


dip 


2 2" 22 
p ft RJ p ft RJ 


(C 1) 


The corresponding sectional lift is 


p ft Rj 


l _ tr 0 - E) 
2„3 " " A 


[ p ft 2 R 2 


nfiP ( f 2 * W f 3 


(C 2) 


The integration of the far field, common part and near field terms (bracketed 
separately on the right hand sid^. of equation (C 1) above) will be discussed for the 
piecewise constant and piecewise quadratic schemes. 

Far Field 

To 0(1 /A), this is the field of a dipole line located along the blade midchord and 
can be expressed as the integral of point dipoles distributed along the line. 


^d ip 


yj R, 


i 


d(s/R 1 ) 


b__i r a 

2,->2 J _ o 2 n 3 / j / t ^ \3 


p ft Rj 


J plTRJ (d/Rj) 
1 




( V R i> / 


d(s/R ] ) 


J p ft 2 R 2 (d/R 2 ) 3 


(C 3) 


2 2 2 2 

where d = + (s - r^) . Using integration by parts, this can also be written as 
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^dip 

2 2 
p ft Rj 


(l-o W 


4A 


(r/R^ 


f ± 

J 3s 


2 2 
\P ft Rj , 


(s - r^ds 
d 


(C 4) 


provided g goes to zero at the root and tip. The azimuth variation of g is approximated 
by a piecewise constant variation over each sub-interval of azimuth, while the spanwise 
variation is represented by a piecewise constant or a piecewise quadratic variation, as 
discussed in Appendix A. 

Piecewise constant representation . - 


P 


far 


<y b /R i> 

4A 


( 1-0 



S j+J /R l d(s/R 1 ) 
s/ Rl WR? 


(l-o 
4A (r/R,) 2 



g; 



S=S. , 

J+l 

S=S. 

J 


The function g. which varies with azimuth will have a constant value over the i th 
azimuth subinteVval that will be denoted by g... The velocity induced by the far field is 
given by 


'Hm. = (i. - . O y y g 

ftR A 4 a 4- j 8i i 



a r ( v )/R ii .. I s )* 1 

3( V R l ) L (r 2 d)/R j J I s. 


(C 5) 


where u = s - r^. The integration will be easier if the derivative can be taken out of the 
integral and this is done as, 


e : 


W 


far 

ft R, 


( 1-0 
4A 


9 (v /R ) ^ ^®ii / 
y -* 0 diy o ,K l ; i j l] e. J , 

3 o ‘ i+ 1 


2 2 2 2 
where now d = + (y^ - y ) + u . Noting that 


(y b - y o> u/R l 
(r 2 d)/R^ 


IS; 


d0 


j+l 


S i 


_y_ _ i r i i i 

2 ~ 2 Ld(d - u) ' d(d + u)J 
r d 
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the integration can now be evaluated using the integrals listed in Appendix B with 
y^/R j = A0 > where A is the inflow ratio. 


0 



(y b - y„> U /Rj 

(r 2 d) /rJ 


d0 


" R° *4 +XI 6 


6: 


6, 


i+1 


Carrying out the differentiation and letting y Q -*.0 , the induced velocity is obtained as 


w 


far 
£2 R. 


= £ £ (Aw. ) p,C * g iS 


i J 


where, for i ^ 1, 


far'.. 


w far>.. 

1J 


p,c * _ (l-Q 


4 A 






6; 


s. 


j+1 


0 : 


i+1 


s. 

1 


(C 6) 


Here, the azimuth index i is set up so that 0 ] = 0 is the collocation point on the blade 
and the trajectory ranges over - °° <6 <0. For i = 1, the trajectory segment has one end 
on the blade and the unit induced velocity is expressed as 


(A w 


far*.. 


p.c. 

ij 


HjlD jf 

4 A *21 



s. 

J 


Piecewise quadratic representation . - For this case, it is more convenient to 
express the field of the dipole line in its second form (eq. (C 4) ). Using the 
representation of Appendix A as 


) = C 2j + 2C 3j (S/R 1* 


P far 

2 2 
p ft Rj 


!LzJ) 

4A 


V R 1 f V/r 

(r/Rj) 2 j=l (i 2i 



uds/R 2 

d/R 1 
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P far (1 -g) 'i 

_2 D 2 = 4A < 

p Rj ] 


|[{ C 2j * C 3j (r b * s) /R 1 } ~T 


- C 3 i <y b / R l> )l n{(ut d)/Rj}J ‘ + 


The induced velocity can then be written as 


0 ■ 2 

T r =-%PZ E (C 2j ♦ C R * Cj. s/R,) / -i- f^l, 

1 4A i j L ] j 3j 11 0V X aVR^LoVR/J 


* 6 .X 


+ (c 3- R ). [ Q 3 |~ y b d/R l 1 

’ v , Lwr/J 


i+1 J b 


6 i 

-3j , i JaTOn [ <y b /R i ) ln { !ut d)/R i} ]de 

1+1 ° 1 


(C 7) 


It may be noted that 


3 / ., 2 , ,2 2 . 4 2 , 2 . 

9yT ( yb d/r > = <x b - y b’ d/r * y b 1 r d 

y b 


i r_L i y b 

- 2 [d+u + d-u - d(d + 


" 2 L d +u d-u d(d + u) 2 ’ d(d . u) 2 J 

Using this relation, and consulting the integrals listed in Appendix B, the induced 
velocity can be expressed as 


%-ZZ <Aw far ) ; • g,(r ,) 
1 I j 1J J 


(C 8) 
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where g (r.) is the midpoint value of the j spanwise segment within the i azimuth 
interval. 1 In terms of the quadratic coefficients, equation (C 7) can be written for i ^ 1 
as 


w 


far (1 - 


PV E [<C 2j * C Rj . C s/R,), 


a " 4a ^ j L 2 i 3 i 


x(lf 6 - + < c 3 j R s ) i (l 17 - x 2 , I 3 ,/2 


(C,;) : If,] 


,f 


'3j i 18 


0 i 


6 i + l 


s. 


M 


s. 

J 


and for i = 1 as 


= ? f <C2 i * C3 ' r ‘ + C 3i s/Rl>1 '‘ 9 


+ (C y I 20 - (C 3 .)j I 21 ] 


s. 


j + l 


The quadratic coefficients can be expressed in terms of the midpoint and end point 
values. In turn, the end point values can be expressed in terms of the midpoint values, 
as shown in Appendix A, to obtain a final expression for the unit induced velocity of 
equation (C 8). The expression is lengthy and will not be reproduced here. 

Common Part 

The common part is the pressure field to which the far field tends at points very 
close to the dipole line, when it behaves essentially as a two-dimensional dipole 
corresponding to the doublet strength at that spanwise location. 


p 

r common 

2 2 

pJTRJ 


(1 -3 y b /R l 

2A (r/Rj) 2 


g(r b ,e) 

2 2 
p JTRj 
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The induced velocity due to the common part is calculated as 


w 

common 
ft R| 


/ 


3 

a(y b /R 1 ) 


common i 

2 2 
p ft RT 


/ 


de 


(C 9) 


Piecewise constant representation . - Noting that 


-JL (y A- 2 ) = ( X 2 _ y 2 ) / r 4 

3y b K b V ' 

equation (C 9) becomes 



The index j represents the spanwise segment within which the i trajectory segment 
lies; in general, the trajectory segment may range over more than one spanwise 
segment, in which case the azimuth interval must be subdivided so that each of the 
subintervals lies entirely within one spanwise segment. 


w. 


common 
ft R, 


= E (A 


p.c. 

) - g- 

common l °ij 


w 


(C 10) 


where, from Appendix B, 


p.c. 

( A w^ ) 

common . 

l 


- [\ 2 '2 * 2 X iV C 3 * (X 5 * X2,, «]| e ' 

i i+1 


Piecewise quadratic representation . - The spanwise variation of g is given by 
equation (A 4) of Appendix A. The induced velocity is given by 
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w 


common 
R , 


(1 - £) 
2A 


i As 


M 


(R. 


Vi ,(R i 


- V c i 


* R s (2R i- r j- Vl )C 2 + R s C 3 } * Si,i + l{ (R i- s i )<R i- r j )C l 

* R s (2R i - r j ■ V c 2 * R 5 c 3}- 2 *i <r j , { (R i- s i )(R rVi )c i 

* R s (2R i- s |- Vl )C 2 + R s C 3 }] 

where 


■[ 


X?lS + 2X. X £ + (x 2 
i 2 i s 3 s 


X?£ + 2X. X ij + (X 2 
l 3 l s 4 s 


X 2 ) ij 


A 2 ) 1^ 


tf. 

]i:. 


Cj = [xflj , 2XjX s l| * (X $ 2 - X 2 ) l|] | e ‘. + 1 


Both g.. and g. . ^ can be expressed in terms of the midpoint 
velocitV can be 1# nnally written as 


values, and the induced 


w 

common 


n Rj 


L 1* 


w >p*^ 


common 




1 ) 


(c 11) 


Near Field 

The near field is a distribution of two-dimensional dipoles along the chord, of 
intensity proportional to the surface pressure differential. 


■near 

2 2 
pfl Rj 


lb 

ttR 


x 

/ 


1 -1 


-^2 4M * ( F 2 + <r b /R 1 )F 3 |C^ 

(x - Xj^Rj/b) 2 + (y b R A /b) 2 


dx 
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where b is the semi-chord. 

To mak e analyti cal in tegra tion possible, the chord is divided into K segments and 
the factors Vl-x/l+x and Vl-x are replaced by an average value over each segment. 
The chordwise integration can then be carried out, resulting in 


8(y b^ R P \p A? 


near 


K f. 
T ± 
k=l 71 


( 2A x " X b) /R l 


n k+ 1 


L{(lA x - x b) 2 + y b} /R i 


where f, 


n 2 o 2 k 
p ft Rj 


HfP [ F 2 + (r b /R l )F 3] b k 


4> - sin 4 ) 
x, 


4> 


k+1 


4>. 


4>/2 - (sin24>)/4 
A x. 


p k+ 1 


cos '•’k * Ax k = Vl - x k 


Piecewise constant representation . - 


w 


near 


ftR, 


u 

-/ 


9(y b /R 1 ) 


r near 
? ? 
p!J RT/ 


d0 


-Z z\ \ 

‘ k ®i+i 


9; /-bi 


("2XRj x ‘ X i) 

d X 


- X 0 
s 


Vl 


d 0 


(C 12) 
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Here, the factor + ^3 3 * s ta * <en t0 a constar| t for each azimuth interval 

i. Using the integrals listed in Appendix A, the induced velocity can be written as 


W npar — , I P* C TT P»c. 


where 


(Aw 1 ) P = — (y ■£- X.)l? - x l5 
near tt L\2A Rj 1 / 1 s 2J 


0 . 


e. 


k+l 


i+1 k 


<-„LC 


(C 13) 


[(^ -*.)*? -X.?] I *' 

1 1+1 


k+l 


F 2 and contain terms proportional to blade motion, namely, 0 , a Q , a^, b^ and the 
twist e • 0 0 

Piecewise quadratic representation . - Here, the quadratic expression for the 
variation of g over a spanwise segment is used and the induced velocity becomes 


w 


near 


ft R 


2a, 

£ £ k 


1 • i k tt A Sj 


[sij{ (R i- 5 M> <R i • r j ) n o 


+ R (2R. - r. - s. ,) n, + R 2 n 0 1 + g. . . i (R. - s.) (R. - r.) n 
s 1 j j+1 1 s 2 ) °i,j+l( x J 1 ) o 

* R s (2R i- r r s j )n i * R s n 2 I 

- 2g i (r j ){(R i - Sj )(R i -s )+1 )n 0 ♦ R S <2R. - Sj - s jtl > n, 

2b 

* R s" 2 }] * ? ? ~h [ F 2 * V R l ,F jl "o 


i k tt A s. 

J 


where 
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{(& ^ - X ,)'I - x,.»] 


0 . 


e. 


i+ 1 


k+ 1 


) I" 

- x l"l 

e , | 

/ 2 

s 3 J 

e i+ i 1 


k+I 


[(m f, - x i) '5 - x s l l] 


e. 


6. 


i+1 


k+1 


As before, the induced velocity expression can be symbolically written as 


w 


near 


ft R 


= L £ 

1 i k 


^w 1 ♦ U U ) P ' q 'l 

_V near/ ik 6i ) \ near / ik J 


(C 14) 


Limiting Behavior 

The behavior of the expressions for the far field and the common part, in the 
vicinity of the lifting line, will be studied for both piecewise representations. 

Piecewise constant representation . - It is sufficient to look at the expressions 
used for the trajectory segment immediately adjacent to the collocation point on the 
blade, viz, the segment with i = 1. Let the collocation point be located at the center of 
spanwise segment m. The unit induced velocity for the far field can be written out in 
full as below, for i = 1 (see eq. (C 6) ). 


(aw far> 


p.c. 


Ij 


As 0-»-O, d -»•(#- 
R 1 





indicating that the first term in the expansion for d will give rise to a singular term, 
which can be separated out. 
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Lim X, (Aw f ) 
0-*-O Lj= 1 lar lj 


(X 2 - x 2 ) J lir- R il S i +1 

(1- £) U s ' 1 v* R i 


_ i v i 

2 i h (r, - r .) 8li 


(1 - 5) <X s - ^ 1 


2 0 8 lm 


since R. = r . Similarly, by looking at the behavior of the common part expression for 
i = 1, it is found that (see eq. (C 10) ) 


Lim ( A w ) 

q > q common | 


. LLlI) , ,y2 ,2,/ _1_\ 

2A <Xs X> (- r 2 e )8 lm 


since = 0. It can be seen that the two terms cancel each other exactly. 

Piecewise quadratic representation . - In a manner similar to that above, the 
singular part of the far field expression (eq. (C 8) ) can be separated out. 


Lim 
0-»- 0 




(X 2 X 2 ) 

(1 - 0 s 


-x 2 > r ( 

~2 ? ^ C 2j + C 3j R i + C 3j S/R l } V 

r 2 J L ' 


V R ‘ 


ft - R j) 

5 l^- R il ln9 


+ C„. R Ik- - R. I Une 
3j s 1 Rj l 1 


Now, C 2 j = g'(r.) - 2(rj/Rj) 
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Using, this relation and carrying out the summation with j, it can be shown that the 
above limiting expression simplifies to 


Lim 
0 + 0 



i= 1 -I sing 


( 1-0 


1A 


— [ 2 8< r j| - 2R s8' (r m )ilne ] 


The limiting value for the common part (eg. (C 1 i) ) can be shown to be 



Since 

S'< r J = <8i, m+ i - i 

it can be seen that the singularities in the far field and common part cancel exactly. 
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APPENDIX D 


EVALUATION OF CHORDWISE AND SPANWISE 
APPROXIMATIONS 

To facilitate the analytical integration of the pressure field in calculating the 
velocity induced on the blade, approximations are used for chordwise and spanwise 
variations, as outlined in Appendices A and C. This Appendix attempts to evaluate the 
accuracy of these approximations. 

Chordwise Approximation 

This consists of replacing the factors Vl-x/l+x and Vl-x^ by their average 
values over each chordwise segment. Since the first factor is the significant portion of 
the chordwise pressure distribution, the following problem is considered. Calculate the 
v elocity in duced by a steady, two-dimensional pressure distribution (proportional to 
Vl-x/l+x ) on a fluid particle travelling parallel to the chord, using both the exact and 
approximate pressure distributions. 

The (x, y) axes are centered at the midchord with the x-axis being parallel to the 
chord. The exact nondimensional pressure field of the airfoil is given by 



where (tij 4>) is an elliptic coordinate system with 

x = cosh n cos <J> 
y = sinh n sin <j> 

If the trajectory considered extends from x = -2tox=2ata constant distance y, the 
nondimensional induced velocity is given by 


v = 



dx 


1 

2 



sinh n cos <(> + sinh n cosh ncos 2 <f> 

2 2 2 
(cosh n - cos 4>) (cosh ri + cos <(>) 


dx 


(D 1) 
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This integral can be evaluated by numerical integration. 

^ Using the approximation, the average value of the pressure distribution over the 
k 1 n chord wise segment is given by 


P _ 1 

k " x k.r x k 


k+1 

/ 


dx 


k+ 


r x k [ 9k+ 


1 °k 


0l, - sin e k+1 + sin Q 


,] 


where x^ = cos 0 


P e (x, y) 


X f - 


2 TT 


k=l k 


7 ‘ 


dx' 


, , v2 2 

(x' - x) + y 


3 Pc 


3 y 


(x' - x) 


f £k 

k=l 2n (x'-x) 2 +y 2 


k+1 


v = 


K p, 

r inn 

kti 4,1 


[ ( x kt i- x) 2 .y 2 i 

( (x k - x) 2 * y 2 ] 


x = 2 


x= -2 


(D 2) 


Two approximations are considered: (1) three segments along the chord, A x^ = 0.1, 0.9, 
1.0; and (2) five segments along the chord, Ax, = 0.1, 0.4, 0.5, 0.5, 0.5. The results for 
the induced velocity are listed in table D 1. Tnese results are also plotted in figure 3 
and it can be seen that the approximate results are quite close to the exact results, 
even down to a vertical distance of 1% of the semi -chord. 
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TABLE D 1. - COMPARATIVE EVALUATION OF CHORDWISE APPROXIMATION 


Vertical 

distance, 

y 

Exact 

Three 

segments 

Five 

segments 

.01 

.5773 

.5708 

.5738 

.05 

.5766 

.5702 

.5731 

.10 

.5745 

.5681 

.5711 

.20 

.5661 

.5602 

.5629 

.30 

.5529 

.5476 

.5501 

.40 

.5357 

.5313 

.5335 

.50 

.5157 

.5121 

.5139 

.60 

.4938 

.4910 

.4925 

.70 

.4709 

.4689 

.4700 

.80 

.4476 

.4464 

.4472 

.90 

.4246 

.4241 

.4245 

1.00 

.4022 

.4022 

.4024 


Spanwise Approximation 

As discussed elsewhere, the span is divided into 3 segments, over each of which 
the variation of the function g is replaced by a piecewise constant or piecewise 
quadratic function. To evaluate this approximation, the following problem is 
considered. Calculate the velocity induced by a finite wing on a fluid particle travelling 
with the freestream flow past the wing, at various vertical and spanwise locations. The 
nondimensional loading along the wing span is taken to be 


g 


e 



with the z axis originating at one tip of the wing so that z = 1 represents the other tip. 
This type of loading has resemblance to the typical loading on a helicopter rotor blade. 

For a particle travelling with constant velocity in the x direction, at constant 
values of y and z, the nondimensional vertical velocity induced by the wing may be 
written as 


v 




/ 



VT 


2 

r dg 

2 i 

+ y + (z - 5) 


172 


dx 


4 Av 







[ 


JiL 


[x 2 + y 2 + (z-,) 2 ] 3 / 2 '[x 2 + y 2 + (z-,) 2 ] 


1-5/2 
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X 


2 


4Av 



d C 


[ 


(y 2 + (z - c) 2 )d 



y 2 d 


2 

¥-2L 


"y 2 + (z - ,) 2 ]V [y 2 + (z-,) 2 ] 2 d 


■] 


|X=X- 


'x=x, 


(D 3) 


where d = £x 2 + y 2 (z - 5 ) 2 J ^ 2 


The integral with respect to £ can be numerically evaluated. The induced velocity can 
also be calculated with a piecewise constant or piecewise quadratic approximation to g, 
using integrals similar to those listed in Appendix B, but the details of these expressions 
will not be written out here. For the calculation, x^ = - 1 and x^ = 1, while the spanwise 
division for the approximations was done with 2 models: (1) three segments along the 
span, with At; - = 0.5, 0.3, 0.2; and, (2) five segments along the span, with A£ • = 

0.2, 0.1, 0.1, 0. J 1. The results of the calculations are shown in tables D 2 and D 3^ where 
the abbreviations p.c. and p.q. represent the piecewise constant and piecewise quadratic 
schemes. These induced velocity comparisons are also plotted in figure 4. 
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APPENDIX E 

FINAL EQUATIONS AND OUTPUT QUANTITIES 


The basic solution is for the spanwise and azimuthwise variation of the dipole 
strength function, g. With the span divided into 3 segments, the value at the midpoint 
of each segment is given the following harmonic representation over the azimuth. 


g< r j) 

2 2 
pffRj 


N h 

G ° + n L [ G |n cos(n V + 




(E 1) 


In the present calculation, the above formula serves only as an interpolation. 

The boundary condition to be satisfied states that the normal velocity induced at a 
point on the blade by the combined pressure field of all the blades is equal to the normal 
velocity due to blade motion. If w, denotes the normal velocity due to blade motion, 
the boundary condition can be symbolically written as 


w. 


(r bo’\o> * .£„,{? ? Uw far>ij p!i 2 R 2 


8i (r i> 


^1 


no. 
blades 


J 


£(aw ). 

. common i 2„2 
i R, 


♦ ? ?[ U »near>iR ^2 * <^near>J } 


(E 2) 


The contributions Aw can be calculated with either the piecewise constant or piecewise 
quadratic spanwise representation, as described in Appendix C. The index j' refers 
the spanwise segment within which the 1 1 trajectory segment is located (if the i l 
trajectory segment ranges over more than one spanwise segment, Aw must be summed 
over these segments as well). 


gj <r i ) 

pR 2 R| 


G. + 



cos (n V k>) + 


G fn sin<n ’ , bi ) ] 


(E 3) 
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where T. • is any typical point of the i t azimuth interval (e.g., either end or the 

midpoint;. 1 It must be noted that the collective pitch, coning angle and the gyclic pitch 

coefficients are contained both in w, and in the near field contribution (Aw 1 ), which 

d near 1 

is proportional to the functions F£ and F y These are defined below 



For presentation as output, the following quantities are calculated: 

(1) azimuthwise variations of total lift per blade, aerodynamic moment about the 
hub and spanwise center of lift location, 

(2) spanwise and azimuthwise variations of sectional lift, pitching moment about 
quarter-chord and center of pressure location, 

(3) chordwise, spanwise and azimuthwise variations of surface pressure 
differential. The expressions for these quantities are presented below. 


Surface pressure differential, 


4P 
2 2 
pa Rj 



g_ 

2 2 
p Q 



(1 -g) 

2A 


(f 2 + 




where x is nondimensional with respect to the semi-chord. 
Sectional lift, 


- - \ U-0 


2 3 
p JTRJ 


2 2 
L p fl Rj 


( 1-0 

4A 



Pitching moment about quarter-chord, 


m 

pfi 2 R^ 


— 3 (1 -5) 3 (F 2 + F ) 
16A 3 2 R i 3 


Center of pressure location (from leading edge), as a fraction of the chord, 



0.25 + 


A m 
(TTD Rj* 


Total lift per blade, 


L 

pIi 2 R* 




n(l - Q 2 j (1 - z 2 ) 
UA 2 ' 2 


a o * [ 2 “ e o< 1 ' e > 


aj(l -Z 2 ) -2 Me] 


+ [2ya o (l - 0 -b 1 (1 - C 2 )] sin ^ + ^b^l - 5)]cos 2\ 
+ [-2y ai (l -5) -P 2 e] sin 2T b } 


(E 6) 


(E 7) 


(E 8) 


(E 9) 


cos T b 


(E 10) 
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Aerodynamic moment about hub, 


M 

p 


' * “ ' 0 J *1 "(*) 


,10^I 2 |0 -£). + 

A 2 I 3 0 


2 2a. o 

ue o ( 1 -r ) } 


+ ( - 4 - 5 - £ 2 ) J COS ¥ b 

+ jya o (l - K 2 ) - ~Y (1 - z?) sin ^ b^l - £ 2 ) j cos 2 ^ 


[-V ajd - C 2 )--^(l + C) j sin 2f b | 


(E 11) 


Spanwise center of lift location, 


[L M_ 

R| " RjL 


(E 12) 


The spanwise integrals of g, required in some of the expressions above, are given in 
Appendix A. 

If the collective pitch, coning angle and cyclic pitch coefficients are considered 
unknown, additional equations are necessary. These equations are given below. 

Piecewise Constant Representation 


- L G°A Sj ♦ <1 - 5> <1 - C 2 > ] a 0 


A C T 


(TT|) B 
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0 


| r j as j * [rk <* - 5) 0 - 5 s ) - a 0 = 

| G ji r i 4s i * [A a-Od-C 2 )]e o -[^ (1-0(1 -5 3 )] a, 


ill u(-* + K + S 3 ) 


- E C?, r j A s j + [& (1 - 5) (1 - 5 2 )] a 0 - [ ^ (1 - 0(1 - S 3 )] b, = 0 


Piecewise Quadratic Representation 

In this case, the equations will also involve the values of g at the ends of spanwise 
segments; however, these can be related to the central values, as described in 
Appendix A. Let this relation be written symbolically as, 


2 2 
pfl K| 


Then the equations are as follows. 


3 g(r.)/ptt 2 R 2 

= £ h.. — L L 

jti ‘i 4 s j 


] r i g° 1 is. 

■ Sn * * h i + i»j ) *t-\ ~ 

* [sX l* - 0 tl - f2) ] a o = (TTp T 

• E H G ° * E <», h.. ♦ s. . h. . .) 5] A A 
i=l L 1 1 jti 1 l i 1+1 1+1,/As.J 6 


[l?A U-5MI -5 3 ) - y J a 0 = 0 
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(b) Piecewise quadratic. 

Piecewise continuous representations of the spanwise 
distribution of doublet strength. 





Nondimensional induced velocity 
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Vertical distance, y/b 

Induced velocity comparison for piecewise constant 
chorwise pressure distribution. 
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.8 


Vertical distance/span 

(a) Three spanwise segments. 

k. - Induced velocity comparison for 
for spanwise approximations. 
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(c) Azimuth angle, 60°. 


(d) Azimuth angle, 90°. 


Figure 8. - Sectional lift versus spanwise location 
for Case 1, /x = 0.29. 
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Figure 10. - Computation with constant azimuth spacing of 5 , 
Case 2, /x = 0.29. 




2. - Computation with constant azimuth spacing of 5 , 
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representations the computational time required for the execution of a typical flight condition has 
been reduced by an order of magnitude with respect to the requirements of the original analysis. 
Airloads which have been computed using the discretized method for a two-bladed model rotor and 
a full-scale four-bladed rotor are in close agreement with measured results and airloads from the 
original asymptotic analysis. For conditions characterized by significant rotor/wake interaction 
the piecewise constant representation requires a reduced azimuth spacing to maintain acceptable 
accuracy. 
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